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Abstract 

We apply the method of the Renormalization Group (GR), following the Polchinski point of 
view, to a model of well developed and isotropic fluid turbulence. The Galilei-invariance is 
preserved and a universality behavior, related to the change of the stochastic stirring force, is 
evident by the numerical results in the inertial region, where a scale-invariant behavior also 
appear. The expected power law of the energy spectrum (g~3) is obtained and the Kolmogorov 
constant Ck agrees with experimental data. 
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1. Introduction 



The incompressible fluid turbulence, at very high Reynolds number (TZ), is described by physical 
systems with a large number of degrees of freedom that are characterized by a substantial scale 
invariance. These peculiarities are investigated in Statistical Hydrodynamics by models in which 
the Navier-Stokes equations are supplemented with stochastic force terms and the analogies be- 
tween scaling at high-Reynolds-number turbulence regime and equilibrium critical phenomena 
have been proposed. Attempts have been made to apply the powerful renormalization- group 
(RG) methods, which have proved so effective and illuminating in critical phenomena, to tur- 
bulence scaling De Dominicis and Martin |^ analysed these models near the critical 
dimension with a perturbative approach. Yakhot and Orszag ^ have extrapolated the pertur- 
bative RG study from the critical to physical dimension and, on this basis they predicted, in a 
long inertial range, various scaling laws, such as the Kolmogorov energy spectrum Q, including 
constant factors as the Kolmogorov constant, the Obukhov-Corrsin constant and the turbulent 
Prandtl number. But the numerical success of the Yakhot and Orszag's theory does not rest on 
firm theoretical grounds. Indeed the critical dimension (d = 7) is very far from the physical one 
{d = 3) and it is not evident that an e expansion is the best procedure |^. Another difficulty of 
the Yakhot and Orszag's method is related to the breaking of the Galilei invariance produced 
by the Wilson-Kadanoff RG analysis in the presence of an infinite number of marginal diagrams 
For these reasons the perturbative RG approach requires some ad hoc hypothesis [|| and 
produces only a phenomenological description of turbulence |5|,[^],0. 

The goal of this paper is that of reconsidering the RG approach for the homogeneous isotropic 
turbulence, for an incompressible fluid, at very high Reynolds number, but in a non perturbative 
context. The point of view is very similar to that introduced by Polchinski ||8|, |Q, [0] in the 
context of the Quantum Field Theory (QFT). The analogy with QFT follows from the description 
of our model in terms of a path-integral representation of the probability generating functionals, 
using a so-called Martin-Siggia- Rosen (MSR) action |11]. 

The dynamics of the model is characterized by the presence of two scales; one of them corre- 
sponds to the macroscopic size L of the system, and the mechanism which maintains the energy 
stationary is operating at this scale. The second, the Kolmogorov scale ry^, is an internal, typi- 
cally smaller scale, at which the energy, after its transfer by a cascade process 0, is dissipated. 



When the dimensionless quantity — = TZi is very large the system admits a wide inertial range, 



3 

= 7? 

where its statistical properties are universal, homogeneous, isotropic and self-similar. As in 
the Yakhot and Orzag's paper IH, the model describes true turbulence only in this long inertial 
range and the manifest scale invariance of this region is interpretable as the existence of a critical 
phase in the infinite volume limit. 
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The analysis of this paper is centered on the study of the flow equations for the effective 
action of an infrared regularized theory. This method involves an infinity of coupled differential 
equations and necessary requires an approximation scheme. The perturbative approximation 
1 12], for the previous arguments, preserves an infinite number of terms and, in our case, no 



summability criterion of the perturbative series is known. In other schemes, as the derivative 



expansion approximation [13|, |14] the correlations functions are strongly dependent on the 
particular form of the stirring force also in the long inertial region. The problem of the approx- 
imation has its solution in the fact that, in the momentum space, the size of the spectrum of 
two point stirring force correlation function is of order O (x)- From this property, using the 
same flow equation, it is possible to see that, in the inertial range (and only in this region), some 
proper correlation functions are smaller then others. 

The paper is organized as follow: 

In Section 2 we summarize the model concerning the Navier-Stokes equation with a stochastic 
force term following the Martin, Siggia and Rose point of view. 

In Section 3 we analyze the Galilei-invariance of the model and we write the Ward Identities 
for the functional generators. 

Section 4 is devoted to the equation flow referred to the scale parameters which characterize 
the model; we also analyze the boundary conditions and some properties of the vertex functions 
are discussed in Appendix A. 

Section 5 contains a discussion concerning the approximation scheme with the relevant hy- 
potheses and a plausibility argument which further considered in Appendix B. 

Some numerical results, concerning the renormalized viscosity and the Kolmogorov constant, 
are contained in Section 6. 

Finally some conclusive remarks are contained in Section 7. 



2. The model 



Let us consider the Navier-Stokes equations, supplemented with a stochastic force term, 

^y{t, x) + n^(t, x)^u^{t, x) = --j^p{t, x) + u^^ u^{t, x) + r (t, x) (2.1) 

Ot OXP p OXa 
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and with the incompressibility constraint 



^ = 0. (2.2) 



dx' 

The fields are the three components of the velocity field and p is the pressure field; the effects 
of the boundary conditions over the zero mean value field are summarized by the stochastic 
source which is a Gaussian random force with mean zero and covariance 

(r(t,f)/^(t',y)) =F"'3(t,t';f,y) =2P-'3(vJiV(x-y)5(t-t')- (2-3) 

-P"^(Vx) is projection onto solenoidal vector fields, required to maintain the incompressibility 
constraint (p.2|). Since the stirring force /" should mimic the instabilities occurring near the 
boundaries, its correlation length must be of the same order of the system size. In other words 
if L is the length corresponding to the size of the system and 

(r q)f{q", q)) = q"; q, q) = 2{27r)^ P''''iq)hiq)6\q + q')Siq' + q''). (2.4) 

is the Fourier transform of (|2.3D , the function h{q) must vanish for g > = m. The explicit 
form of this function is not very important if the self-similarity behaviour is realized, but in the 
stationary limit the energy dissipated by the viscous forces must be compensated by the energy 
introduced through the boundaries and a physical normalization condition is needed, e. g. 

where £ is the rate of power dissipated by a unit mass of fluid and can be calculated by the 



local quantity [0],|15] 



1 2 

£{t,x) = -U {daUj3{t,x) + di3Ua{t,x)) . 

The probability generating functional is given by the functional integral 

W (j, J,p,p) = J PnPuP/e'/'^^^("'™^'-^'^)e~^/'^"'^*^"(")^^^^^ (2.6) 
where we have set x = {t,x). L is the Navier-Stokes density of Lagrangian 

C{x) = n"(x) ( ^ - i^V^ j Ua{x) + u°'{x)up{x)^--^Ua{x) 
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+ ij5(x)-^n"(x) + r{x)u^{x) + u^{x)J^{x). 
p ax" p ax" 

The field if is tlie conjugate variable to the velocity field n", p (the pressure field) and its 
conjugate variable p are the two Lagrangian's multipliers related to the solenoidal constraint 
for the fields and and finally J° and are external sources. In ( p.6[ ) it is possible to 
integrate over the functional gaussian measure 'Dfe 2 J ap-i ami to reduce the space of the 
solutions only to the independent transverse degrees of freedom. We introduce the transverse 
variables v°' = P"^(v)u^, v'^ = P°'^{'\/)u^ and, after the integration over we obtain ^ the 



Martin, Siggia and Rose probability generating functional |11| 



W{J, J) 



DvDve 



i(Siv,v)+J dxJ°'{x)vaix)+J dxv°'{x)Ja(x)) 



(2.7) 



where the action S{v, v) is 



S{v, v) = JdtJ d^xv'^it, x) Udt - u^^) Va{t, x) + v^{t, x)dpv''{t, x) 
+ ydtjd^xj d^yv'^it, x)N^fs{x - y)v^{t, y). 



(2. 



The multitime and multipoint correlation functions are obtained by functional derivatives of 
(|2.7|) with respect to the external sources J" and J". The connected functional is defined as 



lnW{J, J) = iZ{J, J) 



(2.9) 



and the generic connected correlation function is given by 

(xi)...^;"" (x„)?)""+^ (x„+i)...i)°"+"' (x„+™)). 
5 S 5 S 



n+m —iZ{J,J) _ 



(-i)"+"e 



6>^{xi)"'5>" {Xn) SJ^^^+^iXn+l)'" SJ''"+"^{Xn+m) 

where the stability conditions 



JZ{J,J) 



(2.10) 



J=J=0 



6Z{J, J) 



(5J"(x) 



6Z{J, J) 



J=J=0 



(5J"(x) 







(2.11) 



J=J=0 



are satisfied. Finally the convex vertex functional (the effective action) T{v,v) is defined by the 
Legendre transformation Q 



Z{jJ) = r{v,v)+ I dx(j"(x)wa(x) + i)"(x)J«(x)) . 



(2.12) 
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Since from ( |2.12| ) we have the relation 




5T{v, v) 




the stabihty conditions ( |2.11| ) are rewritten 



6r{v,v) 



u=u=0 



ST{v, v) 



u=u=0 



= 0. 



(2.13) 



We conclude this section observing that the knowledge of the functional W{J,J), and then of 
Z(J,J) and T{v,v), gives a complete solution of the model in terms of correlation functions. 
The direct calculation of the functional generator from the Navier-Stokes equation (|2.lD , with 



suitable asymptotic boundary conditions, is a formidable task which may be simplified studying 
the constraints deriving from the symmetries of the model and the scaling properties connected 
to the variation of the dimensional parameters. The following sections are devoted to developing 
this point of view. 



3. Galilei invariance and Ward's identities 

In the theory of turbulence we consider the very important hypothesis that, in the limit of 
infinite Reynolds number, all the symmetries of the Navier-Stokes equation, that are broken by 
the mechanism producing the turbulent flow, are restored in a statistical sense at small scale 
and away from the boundaries [^]. In our model, the last term in eq.( p.8D is Galilei invariant, 
as the original Navier-Stokes equation, and this symmetry will be transferred to all correlation 
functions. Then, as a physical system is described by its symmetries, if scale invariance is also 
recovered in the inertial range region, the model will provide a good description of the isotropic 
turbulence in this region. In this section we study the constraint deriving from Galilei-invar iance. 

Let us consider the following coordinates transformations 



.OL 




with corresponding field transformations 



y'°(t,f') = y"(t,f)-c' 
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where c" is a constant velocity. We have then the Lie's derivatives 

5v"it,x) = -cHdxv"it,x) - c", 
6v°'{t,x) = -cHdxv°'{t,x). 

Prom (|3.1| ) we obtain the following Ward identity for the functional generator W{J,J) 



(3.1) 



g^wijj) = I dx \it,r[x)S^S— + itr{x) ^ ^ 



1 >V( J, J) = 0. (3.2) 



5x^5J"(x) dx^5J'^{x) 



The functional-differential operator Q'^ carries a representation of the infinitesimal Galilei trans- 
formations and it is a generator of a closed Lie algebra. Indeed the model is also invariant under 
time and space translations, given by the Ward identities 



nWiJ, J)= I dx (iJ"(x)|- / ^ + iJ°(x)-^ ^ 



dt5J^{x) 'dtSJ<^{x) 



>V( J, J) = 0, (3.3) 



V'W{J, J)^ Idx iJ-{x)^-^ + iJ-(i)A^^—] W(J, J) = (3.4) 



dx^ 5J^{a 



which obey the commutation relations 



g^^n =-v^, g^,g^' = g^,v^ = v^,v^ = v^,n =o. 



(3.5) 



Analogous relations hold for the functionals Z(J,J) and T(v,v). Introducing the Fourier trans- 
form 



1 



(2vr) 



2 , / (i3ge-^«°*-^«^^\"(x),... 



(3.6) 



the previous relations are rewritten, for example in terms of the vertex functional, as 



(3.7) 

-nr{v,v) = |d^|^gV(g)^^ + <7°r(g)^^^r(z;,t)) = 0, (3.8) 
-r^riv,v)= fdq(q^r{q)-47^ + q'r{q)-^]T{v,v)=0. (3.9) 



Sv'^iq) 



6v (q)^ 
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In order to derive from the previous identities, some general properties of correlation functions, 
we introduce the following simplified notation for the n-points vertex functions. 



5"+'"r(w,{;) 



(5i}°i {qi)..6v°'" {qn)dd°'"^' (g„+i)..(5t}°"+- (0) 



v^l {gi)..v°'" iq„)i, "+i {q„+'^)..v°'r^+m (q) 



(3.10) 



v,v=0 



We note that the last field in ( 3.1C| ) is written in configuration space, consequently only the 
independent momenta appear and no kinematic 5 singularity is present. 

The first non trivial identity is 



0, 



(3.11) 



and a consequence of this relation is, for example, that the coefficients of the two terms of the 
convective derivative dt + v^da have the same value. Indeed from power counting and tensorial 
analysis we set 



and 



r..(,)..(fc)..(o) = ^ (^^V + 5^"^') (^7^ + S(p, k)) 
with n(6) = S(6,6) = 0. From (|ll|) we obtain the constraints 



d 

o-Q = (^TZ and -^U{k) = iS(0, k). 



(3.12) 



Analogous relations follow from the general Ward identity of the proper vertex correlation func- 
tion that can be written as 



k'^T 



A ^ 

'q^'- v°'{k)vi{p)...v>^{0) '^P 'Q^'^v°'{k)vt{p)...vl^{0) + ••• + ■'^{j^(0)C°'(fc)CT{p)...i)'3(0) ~ (3.13) 



4. The evolution equation 



The mechanism of energy input in the system, necessary to preserve the stationary regime, 
is described by the last term of the Martin, Siggia and Rose action (|2.8|) and it contains a 
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characteristic physical scale. As we have already observed this scale is given by the parameter 
m that is the width, in momentum space, of the two points stirring force correlation function 
support and m = where L is of the order of the linear dimension of the system. In a real 
physical system an other short-distance scale, related to the dissipative region, is present; the 
relevant parameter is the Kolmogorov's cutoff = — . The two scales m and Kd are very 
important in the physical description of the system. In fact the dimensionless quantity 



m 



n (4.1) 



is the Reynolds number |^ which is a measure of the nonlinearity occurring at the scale L. 
If the Reynolds number is very high the system develops a wide domain of the Fourier space 
{the inertial range) in which the non linear terms dominate. If we consider the momenta in 
a region included in this domain, such that m <C g <C K^, both boundary and viscous terms 
can be neglected and we expect the system to be strongly self similar, homogeneous in space 
and isotropic. This behavior of the system is interpretable as the presence of a critical phase 
in the limit IZ ^ oo. In order to reproduce the previous phenomenology starting from the 
model described by the action (^.g|), we consider the Kolmogorov cutoff as a reference scale 
and we study the system assuming the other scale m to be very very small compared with K^- 
In other words we study the system near the infrared fixed point (if this fixed point exists) 
connected to the previous described critical phase. We expect that, for momenta q such that 
m <ti q <^ Kd, the self similar behaviour becomes manifest and that correlation functions obey, 
with good approximation, the behaviour given by dimensional analysis. In particular the energy 
spectrum E{q) must follow a power law of the form Q 

E{q)=CK£h-~^, (4.2) 

where Ck is a constant of the order of unity and the renormalization group methods allow a 
quantitative estimation of the Ck Kolmogorov constant. Experimental tests prove that the 
scaling behavior of correlation functions of order higher then 3 has some corrections. Such 



corrections are generally linked to the intermittent behavior of turbulent systems [16|, i. e. to 
strong, nonlinear, rare events that are usually associated with the instantonic configurations. 
This aspect is not discussed in the approach here analysed. 



In order to obtain an evolution equation for the model (2^) we introduce the six components 
field ^"■{q) = (y°' (q) , v°' [q]^ and the density of lagrangian 

£((?; m) = C\q) + U>'^{-q)R,,{q- m)^>\q) (4.3) 
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where J~P{q) is the original Navier-Stokes lagrangian (independent of the scale m) without stirring 
force, which is totally included in the last term where the matrix RabiQ] n^) is given by 



Rabiq;rn) 





2ih{q;m)P°'l^{q) 



(4.4) 



Concerning the function h{q; m), it must satisfy the physical constraint (|2.5| ) and its momentum 
spectrum should be of width m. A possible choice is 



h{q;m) = f— ) , 



(4.5) 



where, for example 



X 



m 



^ ( — ) e rn^ 



Dq is a dimensional parameter fixed by the constraint (2.5), once the homogeneous function 
xim) specified. The function ( |4.5| ) is very similar to that chosen by P. Brax [17| and as the 
parameter m is different from zero, the density of lagrangian (^^) corresponds to a model with 
a Galilei-invariant infrared regularization. The ultraviolet convergence is assured by the same 
function ( [4.5D and by the existence of a physical cutoff Aq > where, for momenta g ~ Aq, all 
statistical fluctuations rapidly vanish due to the dissipative effects. 

Taking into account that all the m dependence of the functional r(^'; m) (r(^'; m) = T{v, v)) 



is contained in the function h{q;m), one easily M, |12], |ffq] derives an evolution equation in m 



mdm { r(^;m) - ^ / dq^f''i-q)Rab{q;m)^''iq) 



dq 5'^T(^;m) 



(4.6) 



where A^,^^^(g;m) is the full propagator which, in term of the two points proper correlation 
functions, is given by the matrix 



V 



P°"^(g) 





(4.7) 



The auxiliary functional 



(5^r(^;m) 



Tab{q, q'] m) is defined by the integral equation 



fab{qA';m) = T2\q,q';m) - / dg'T ,,(-g", g'; m)A^,*,(g"; m)C*(g, g"; m 



(4.8) 
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where T^Jl^{q,q';m) is obtained by the functional 

r-* i^; m)=T{^;m)-^J (i#,(-g)r^^(^)^^(o)§fe(-g), (4.9) 

i, e. the generator r(^';m) without the bilinear terms. In eg. ([4.7]) we have assumed that 
T^ay0 = 0; in perturbation theory all vertex functions containing only field vanish, but 
this property is an exact consequence of ( [4. 61 ) and the boundary conditions (see appendix A). 



From (4^) we obtain the equation for the vertex correlation functions. In order to tackle 
specific problems, consider the flow equations for the two points vertex functions r^a(p)^/3(o) and 



'n)v0(nV^ they are 



m(9mr{)a(p)£,/3(o) — J -j—^mdmh{q;m) (^A£,7„p(g; m)rjp(^)^5(_^){)(.(p){)/5(o) A^a^T (g; m) 
-Ai)7„p(g; m)ri)p(^){)a(j3),^s(o)A,fs^r(g + p;m)r^r(.^^^-a(_^j.^(QjA^a£,7(g; m) 



Ai)7„P {q\ m)r-p( (j3_ A^.^r {p-q- m)T ^ -^^s A^s^^ {q; m) j (4. 10) 



and 



(r-(^),,(o) - 2ih{p-m)P^^{p)] 
-^md^mKq;m) (A{)7i;p(g; m)r-p(.)-,(_.).<.(^).^(Q)A^i^,7(g; m) 



(2vr)4 

-A^7^p(g;m)r-p(^).«(^)^,(Q)A^.^r(g+p;m)r^,(-^^)-,(_.^-^ . (4.11) 

In the r. h. s. of the previous equations there appear three and four points vertex functions that 
are defined by flow equations involving higher correlations thus producing an infinite system. 
We will analyze the fundamental point of truncating this system in another section, studying 
here only the general properties and the boundary conditions. 

If there is an infrared fixed point in the limit m — > 0, it is convenient to study the system 
near this point. For this reason we measure all momenta in terms of the spectral parameter m 
and set 

p^ = Tm?' , \p\=p = xm, q^ = rjw? ^ \q\ = q = ym, ... (4-12) 

We suppose that the spectral parameter m is very small compared with the Kolmogorov cutoff 
Kd. The space of the x-like variables is then divided into three regions. In the first region (the 
infrared region) the p momenta are smaller then m and 

X <l. 
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Here all correlation functions are strongly dependent on the form of the stirring force. In the 
second region (the inertial region) the p momenta are finite and m is very small compared with 
these momenta. Here we have 

X > 1. 

In this region we expect a strong self-similarity behaviour of the system, provided that 



X < 



m 



The last region is the dissipative region, where 



X > 



m 



The boundary conditions for eq.s( 4.10| ) and ( [4.11 ), and for all possible flow equation, are obtained 
in this region by the consideration that for 

X = xo = — ~ oo, where Aq > 
m 

all statistical fluctuations vanish as a consequence of the dissipative effects. In other words for 
any external variable x of the order ^ ~ cxd all correlation vertex functions with a field content 
different from the vertices of the original Navier-Stokes action, without the stirring force term, 
vanish. Therefore the boundary value of the effective action is the classical Navier-Stokes action. 

From equations ( 4.10| ), (4.11) and the previous considerations it is possible to extract some 
general informations about the vertex correlation functions. Let us consider the two point vertex 



functions F,- 



C"(p){)'3(0) 



, that we rewrite, in terms of the rescaled variables and taking into account 



the power counting analysis, as 



-U\r,^]r + f(r,x,^\x' 



(4.13) 



-1 -pm 



2i- 



x{x)+M\ t,x,^ 



P"^(p). 



(4.14) 



These functions are defined, at m 7^ 0, for all values of the variable r and x. In particular they 
satisfy T™ ^.^ o)i)/3(o) ~ ^ M(r, 0, ..)= for any m 7^ 0. It is possible also to see that all vertex 
functions (with the exclusion of T{,{,) are almost proportional to the momenta of the v fields. If 
we go back to equation ( 4.11| ) we see that some inconsistencies are present: indeed the r. h. s. of 
this equation, for x = is different from zero due to infrared divergences of the full propagator. 
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This apparent inconsistency is caused by the fact that in the function h{x; m) 



A) 



xix) the 



value X = is not accessible due to the finite size of the system and therefore we must define 

X(0) = 0. (4.15) 



The equations ( 4.11| ), ( [4.10I ) are defined for arbitrary values of the variable x and the previous 



condition for the function h{x; m) at x = is a physical regularization connected with the 
infinite volume limit. Taking into account condition ( [4.15| ) we obtain a flow equation for the 
quantity 



- {i" (T,0)t)'5 (0) 

















2Tdr) 















{mdm — 2Tdr) I T, 



which has the general integral 



1 = 1\Dq ''m^T 1 . 



But the asymptotic boundary conditions in x = (and t = ^) enforce I I^Dq '^m t J = 1, 
independently of the values of the m variable. Therefore we have the exact result 

/(t,..)^1. 

Thus the two points vertex function is rewritten as 



■pm 



-zr + / I T,x, ^ I 



(4.16) 



From the Ward identity ( p.llD and from ( |3.12D we have another exact result concerning the 
three points vertex function r^A(g)j;Q(p_g){)/3(o) which is 



-p -pm, 

= imx 



(4.17) 



where p°' = mxn°'. These exact results, which are in agreement with the non renormalizability 
theorems of this model, were already known in perturbation theory. 
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5. The approximation scheme 



In order to study the structure functions of the two points correlations functions, let us consider 
eq.s ( [4.1C1| ), (4.11). As we have already observed they are equivalent to an infinite system of 



flow equations; therefore a truncation procedure is needed. The important observation is that 

the correlation vertex functions, which have a different field content with respect to the original 

classical Navier-Stokes action, are relevant in the infrared region, i,e. when the x-like variables 

are near unity, but they are less relevant in the inertial region where x ^ 1. The physical reason 

for this expected behavior is due to the statistical fluctuations induced by the non local vertex 

2ih{q]m)P"^ {q)v°^{—q)v^ {q) which has a very concentrated Fourier spectrum in the long wave 

length region. The non linear terms transfer, by the cascade mechanism, these fluctuations to all 

regions of the Fourier spectrum, but the speed of this transfer depends on the coupling of the non 

linear term with the degrees of freedom of the system; in other words on the Reynolds number. If 

the Reynolds number is very large the statistical fluctuations are swiftly transferred to a spectral 

region where the dissipative effects are dominant and, in order to preserve the stationary regime, 

the corresponding amplitudes are decreasing for increasing momenta. A more specific analysis 

follows from the connection between the parameters Dq, m, the quantities related to the regime 

of the physical system such as the Reynolds number, and the rate of power dissipated by a mass 

unit S of fluid. From the normalization condition (2.5) we have for the dimensionless quantity 
Do 

rrf* 

Do vr2 {£) 



--4 



The Kolmogorov's cutoff is defined by = with u the kinematic viscosity. Recalling then 



(|4.l| ) the following relation is obtained 

D() TT^ .,„n 

m"^ J dyy^xiy) 

As we see the statistical fluctuations vanish for 7?. — > according to the fact that, in this case, 
no transfer of energy in a spectral region different from m is possible and the energy pumped 
into the system at this scale is immediately dissipated (indeed for TZ ^ m ^ K^). 

If we analyze in more detail eq.s ( 4.1C| ), (4-11) we find in the r.h.s. the three and four points 
proper vertex correlation functions 



^vP{q)v6{-q)v°'{p)v0{O)^ ^ vP{q)v^ {-qYv" {p)vf' {0) (^•^) 
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and the two points proper correlation functions appear in the fuh propagators as results from 
(|4.7]). Concerning these last correlation functions we have a first observation about the function 



From ( 4.16 ) we rewrite the function /(r, x, ...) isolating the terms which depend on the 
T variable i,e. 

/ (^^ X- ^) = / (x; ^) + fr (t, X- ^° 



Prom the Ward identity ( |3l^ ) and ( p^D we obtain 



^/.(r,x;^)x^=.E(r,0,x,0;^) (5.4) 

i. e. fr (^T, x; has the same scaling behavior of the statistical corrections to the three 

points classical vertex. Therefore, for the previous argument, we expect that decreases faster 
then /. In a completely similar way we see that all r dependent corrections of the vertex 
functions with n v and m v fields, are connected, by the Ward identities, to the functions with 
n + 1 V and m v fields. Writing 

^) = -^^9 {x; ^) + T^g^2) {x; ^) + ... + r^^H (x; ^) + (5.5) 

from the last observations we have that g (^x; is related to three point vertex statistical 
corrections, g(^2) (^x; is related to the four point vertex and so on. Taking into account the 
previous considerations we can write the flow equation for f{x;-^), where we analyze only 
the leading contributions appearing in the triple vertex. In other words considering the vertex 



functions that appear in the r. h. s. of equation (4.10), i. e. the three and four points function 

^THv,y)v'^i-r,,x-y)vHo)^ ^Tpir,,y)vH-v,-y)v'^io,x)vHoV Separate the leading contribution from the 
statistical correction by writing the vertex F™^ as 

rZ, = + (AF)„„, (5.6) 

and we suppose, for x S> 1, the following approximations 

limF-,~F2 and lim F™,, ~ 0. (5.7) 



The same analysis is performed for eq.(4.11) where the vertices involved are F" 



F™, ,-p, F!!5;, , with the approximations 

v^{r],y)v {-ri,x-y)vP{0)^ v^{ri,y)vP{-'ri,-y)v (0,x)vp{0)^ 



limF:;i^~0 and lim F^;;^^ ~ 0. (5.8) 
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In this way eq.s ( 4.10| ), ( [4.11 ) give a closed system of differential equations. The approxima- 
tion here proposed can also be obtained by a suitable field truncation, preserving the Galilei- 
invariance, in the effective irreducible action ||n 



In order to write, from ( [4.10| ) and ( ^.111 ), the equations for the structure functions / 



and M [ x 



Do 



we define 



^1 



M\x-^]=M^{x) 



and after straightforward calculations we obtain 

(mdm - xdx - (l)m{x) 







(3 


^ydv)x{y) 








(x^+y2+2txy)+<t>miy)y'^ 



X 



. Do 



(5.9) 
(5.10) 



(5.11) 



{mdm - xdx - 3) Mraix) 

(l-t2)a;2 (^l+2t2 + ^+3t|) {3+ydy)x{y) (xm (^^/ x'^ +y2 +2txy'^ +Mm (^^/^^W+2t^'^ ) 



4>m{y)'Pn 



x'^+y'^+2txy ) {x'^ +2txyY 



4>m ( ^/x2+y2+2txy ) {x^+y^+'2txy)+<l,m{y)y^ 



(5.12) 



Concerning the boundary and initial conditions for ( ^.11 ) and ( |5.12| ) we must consider that 
these equations are valid only in the inertial region, i.e. for x ^ 1, provided x <C ^ = Xd- 
As the parameter m is very small compared with Kd (we are studying the system near the 
infinite volume limit) it follows that ~ oo is a good approximation. Recalling that the 
Kolmogorov scale is related to the viscosity and to the rate of energy dissipation by the 

relation Kd ~ (p')* 0]' previous approximation, for m 7^ 0, is equivalent to ~ 0. In 
conclusion we extrapolate the following boundary conditions: 



(00) = Mjnioo) = 0, for m / 0. 



(5.13) 



It is clear that ( |5.1l| ) and ( |5.12| ) with ( |5.13| ) enable us to describe the physical system in the 
dissipative region. Considering that we are studying the system near the supposed fixed point in 
the infinite volume limit, the initial conditions are given for m — > 0. As we have already observed 
our model corresponds to a regularized theory and, for vanishing regularization parameter m, 
the infrared divergences, due to the infinite volume limit, will appear again. The choices (B.9|), 
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( 5.10 ) fix these divergences according to the power counting and dimensional analysis; in (pmix) 
and Mm{x) no explicit dependence on the parameter m is present. It is also evident that for 
X ~ 1 the considered equations are inadequate for the description of the physical system. 

The meaning of our approximation scheme is evident by observing eq.s( p.li ), (5.12) which 
correspond to the first step of the loop expansion. However these equations contain a nonper- 
turbative contribution concerning the structure functions { f).9\ ) and ( 5.1C| ) which are solutions 
of the same equations with the initial conditions for m ^ 0. 

The next step requires the one loop calculation of the vertex functions (Ar)^^^-, T^^^, T^- 
and - - still using the ERG containing in the r. h. s. only the leading terms and the structure 
functions ( |5.9| ) and ( ^.lOl) . The two loops contributions follow inserting these vertices in the 
original equations ( 4.10| ), (4.11). This procedure is extendible to an arbitrary order and it 
contains a non perturbative input given by the structure functions ( |5.9| ) and ( 5.10| ). Naturally 
this scheme is meaningful only if the two loop are smaller then the one loop contributions, and 
so on. We also observe that, in our approach, no particular hypothesis concerning the effective 
action is needed and this result follows from a self-consistency argument of our model. 



6. Numerical results 

In this section we report some numerical results concerning the solutions of the system ( |5.11| ), 



(5.12). The analysis takes into account different choices of the stirring force, all satisfying the 
spectral condition discussed in the previous sections. By the expected universality in the inertial 
range we must verify that the exact shape of the noise correlator is not very important. We 
consider a x(^) function with the general form 

X{x) = x^e-^'" (6.1) 



where s and n are positive integers. Eq.s ( 5.11| ), ( |5.12| ) are solved with an explicit finite-difference 



scheme starting at the boundary conditions (5.13) and by the initial conditions which require 



that no explicit m-dependence is present for m ^ (i. e. for the system near the critical 
point). The numerical computations show that the solutions maintain explicit independence on 
the scale-parameter m. The function cp^x) and M{x) are computed in the range 10~^ < x < 10^. 
The inertial range corresponds to x ^ 1 and in this region the function 0(x) and M(x) reach a 
scale invariant regime. Indeed, referring to the parametrization 



ax-t+" and M(x) = 7x"^+^, (6.2) 



17 



Functions M, <^ 




10"- 



10 



100 



y=p/m 



Figure 1: (j) and M functions 

the numerical computation gives for a and /3 a value ~ +0.37 in the region x > 1 and the 
function (j){x) and M{x) are graphically represented in Fig.(|l|). The detailed numerical values of 
the parameters related to the functions <^(x) and M(x) obtained with different noise functions 
X are summarized in table (1) and the substantial independence of the exponents — | + a and 
— 3 + /3 on the particular form of xi^) appears very clearly. 





cr (±0.001) 


-| + a(±0.0005) 


7 (±0.001) 


-3 ± /? (±0.0005) 




x'^e~^^ 


0.899 


-0.924 


20.00 


-2.590 


14.85 


1 — 

X e 


0.900 


-0.935 


22.49 


-2.601 


43.55 


1 — 

X e 


0.935 


-0.930 


26.02 


-2.630 


52.46 


4 — 

X e 


0.924 


-0.933 


29.63 


-2.599 


5.94 


4 — 

X e 


0.994 


-0.929 


28.97 


-2.595 


42.95 


x^e~^ 


1.02 


-0.930 


33.55 


-2.596 


66.81 




1.03 


-0.930 


35.93 


-2.633 


88.82 




1.08 


-0.929 


38.05 


-2.595 


109.97 




1.09 


-0.931 


38.41 


-2.597 


109.97 


1 2 

e ^ 


0.942 


-0.922 


27.70 


-2.588 


54.86 



Table 1: Parameter of the functions 0(x) and M(x) for various stirring forces. 
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A significative check of our analysis is given by the computation of the energy spectrum. From 



E{x; m) 
1 



1 



6^Z 



2{2ttY 



dq^q^dqdn 



2 5J'^{x)6Ja{x) 
1 



1 



{u°'{t,x)Uait,x))c 



J=0 



dqE{q; m) 



(6.3) 



we obtain 



E{q;m) = ^q 



1 „ M(^) 52 1/7 



1 



^1 _5 
t 3g 3 



2 Vo-/ \! dyy'^xiy) 

Recahing ([4. 21) the Kolmogorov constant Ck is given by the g-independent part of (|6.4D i.e. 



(6.4) 



1 



2 Vo-/ Kjdyy'^xiy) 



(6.5) 



The energy power behavior in terms of the q- variable is given by the exponent ^ = — 1 + — a. 
The numerical values of Ck and ^, computed with different noise functions x(^) are summarized 
in the following table 





C (±0.001) 


Ci^(±0.002) 


X{x) 


C (±0.001) 


Ca-(±0.002) 


x^e~ 




-1.666 


1.124 


X e 


-1.666 


1.461 


x^e~ 


-x^ 


-1.667 


1.267 


x^e~^ 


-1.666 


1.660 






-1.666 


1.417 




-1.666 


1.767 


_ 1 
e ^ 




-1.667 


1.489 




-1.666 


1.784 


4 - 
X e 




-1.667 


1.624 




-1.666 


1.785 



Table 2: Slope and Kolmogorov constant Ck for various stirring forces 

Looking at table (2) the exponent ^ in ( |6.4D appears to be practically independent on the 

5 

3- 



particular form of xi^) and corresponds to the Kolmogorov's value — t. The constant Ck is 



stable for very narrow x functions and approaches the value Ck ~ 1-78 compatible with the 
experimental results These results are summarized, for the energy spectrum E{q;m), in 

Fig.(|). 

We conclude this section giving a numerical justification of the loop expansion here consid- 
ered. A rather rough evaluation in Appendix B shows that the two and one loop terms 



^^-^mam/i(g;m)Ai7^p(^;m)r^p(^)^5i(^^)^5<.(j5)ii/3(o)A„i5^(g;m), 



(6.6) 
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Energy Spectra 




1 10 100 

y=p/m 



Figure 2: Energy spectra for various x(x) functions 

.^^mdmh{q; m)Ai^^p{q; m)TfJ^.^^^^.^.^^^^Acs^r{q + p)r(°^(.^^^_,(_^^.^(Q^ A^,^^(^; m) (6.7) 

are in the ratio ^ i.e. proportional to the spectrum size of the noise function h{x; m). Analogous 
results are obtained for the other two loops terms and in general the terms of our loops ordering, 
generated by the ERG, are decreasing as inverse powers of n. This result is not a sufficient 
argument for the validity of the considered loops expansion, but only a proof of consistency of 
the theory and a precise numerical estimation of all interesting quantities is very hard. We also 
observe that, unlike a naive perturbative loops expansion, the terms of the expansion are not 
increasing as a power of the Reynolds number Ti. 



7. Conclusions 

The numerical results in the previous section are encouraging and, as reported in appendix B, 
an essential ingredient appears to be the size of the spectra of two point stirring force cor- 
relation function. If the approximation scheme here proposed is physically meaningful, the 
non-perturbative contributions are mainly due, in the inertial region, to the two point function 
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^'uu ^-iid, particularly, to the anomalous dimension of the u field (no anomalous dimension is 
associated to the u field as results by the Ward Identity). In other words, near a critical point, 
the ERG provides a summability criterion concerning the two point function and the residual 
non classical dynamics is described by perturbative corrections. We cannot provide a proof for 
this scenario, but only support it with plausibility arguments since a direct check requires a very 
heavy numerical analysis. We think that an interesting alternative to the numerical analysis is 



given by the method recently proposed by T. R. Morris |13], in the context of quantum field 
theory. The Morris's method is similar to our analysis in the treatment of the field anomalous 
dimensions in the ERG, but our perturbative loops expansion is replaced by a more promis- 
ing nonperturbative approximation of the effective action. The Morris's scheme correspond to 
a quasi-local approximation of the effective action obtained by a derivative expansion. How- 
ever the application of this method to our model is not straightforward, in particular Morris 
approach does not preserve Galilei invariance. It should however be interesting to compare nu- 
merical results obtained with different methods in the effort towards a better understanding of 
the different approximations in the description of the physics of these systems. 
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Appendix A 

In this appendix we prove that the vertex correlation functions containing only f " fields neces- 
sarily vanish for an arbitrary value of the parameter m 7^ if these function vanish for a given 
m = m. Consider the flow equation for any such function; from (|4.6D we have 



lll'Um^ i}<^l(pi)...D°'n(0) 

j -^rndmRab{q;'rn)l\,^^^^^{q]m) ^"+"r(<l/;m) 



5^c{q)S^i{-q)&v°'i(pi)...&v°'n.{0) 



*=0 



(Al) 



It is straightforward to see that all terms in the r. h. s. of eq.( [A.lD contain at least one vertex 
with only the v°' field or the component A„a^/5(g; m) of the full propagator. For this reason 



(0) 
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and at the generic point m + 6m near m we have 

rr^t£)....^n(o) = O {i6mf) (A.2) 

which iteratively yields the desired result. In other words we have a non-renormalization property 
for these vertices and they will be absent if they are not present in the original classical action. 



Appendix B 

Here we propose a heuristic argument in order to compare two loops with one loop terms like 
the ones considered at the end of section 6. Taking into account the results of section 6, consider 
the two point functions 

^v^ipmo) - (,,,),.(o) - rn' (-ir + (^) ' x^^^ j P^^in,). (B.l) 

and 

where, from Table (1), oj = 1 + ^— | + ~ 0.07 and f3 ~ 0.4. The calculation of the vertex 
function ^yP(^q)yS(^_^)ya(^p'j{jfj(^Q'^ is performed by the ERG where, as in the case of the two point 
function, only the leading vertex is taken into account. In this approximation the ERG which 
we consider is 

mdm^f^p {p){)/3 (0) 



(0). ^ A ^-nik-v- a-m)T^^^^ ■ 



,-p(0) ^^.„('^_"i \p(o) ^ 



A,,af,d(k — p:m)r^!^].f , + 

" " ^ ^' ^ iX<l{k-p)u"{p)u'^{-k) 



^u\k)u0{-p)u''{k-p)^^''^'^^ P'™'^ 



■^(0) _ A m'nik — v — a"rn)V^^^ ~ 

u'-{k-p)u^(-q)u"^{-k+p+q) ^ f I ^n(^_^_^)„a(^)„<i(_^^g) 



K^uAk - q; "^)ri°4_^)„.(,>c(_fe) + ^%-p-q)uP(q)umi-k+p)^---<^ -p + q-m) 
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+ 



u"{k-p+q)u"{p)u"-(~k-q) " " ^ i-i J u'i(k+q)u^{~-q)u<=(-k) ' ' 
ii"(fc-(j-p)M«(p)M°{-fc+<j) ^ ^' u<i{k-q)u 

We rescale the variables (0,p"), {q^,q'^) and {k^,k^) as 



r^"-" A .^('A' -I- f9- TTi'ir'- -I- 



— g — p; mj 
+ ^A„.^/(fc;m). (B.3) 



myn 



(B.4) 



where n", and n are unitary vectors along the directions of p", q'' and k . We also observe 



that the value of 



is constrained, by the support of the function h{k; m) (see (^) and 
( |6.1| )), into a narrow region centered around ~ m, so that, from ( [B.4[) , y ~ 1. Moreover, 
in order to evaluate the two loops term ( |6.6| ), the variable z has also the same constraint. 
From the last observation it follows that, in the inertial region (where x ^> 1), we can consider 
^ ~ ^ <C 1. Taking into account only the leading terms we obtain from eq. (|B.3| ), after tedious 
but straightforward calculations 



{mdm - Xdx - Zdz - 2z^d^o) ri)p(q)^)5(_^){;c(p){i/3 



(0) 



1 l+cj-p ( 2+s-uj \ (Do\ 



a2 (z")2+a2(-2|)^x2+2'- 



+ 



+ ldyx(.y)0{l^). (B.5) 



Where we have considered xiu) = y^^~^ as in ( |6.1D . A solution of eq.( |B.5D is not obtainable 
with elementary procedures and requires numerical methods, but it is still possible to extract 
sufficient information concerning our purpose with very simple arguments. We simplify (|B.5| ) 
rescaling the variable z^ as 



z' = a 



m 



(B.6) 



and setting 



vP{q)v^{-q)v'='{p)vl^{0) " g^2 



1 l+UJ^ /2 + s-w 



n 



2n 



Do 



3 X 



-(l+3a;) 



G(m,2;,u;°)nJn^P"'3(np) + ... 

(B.7) 
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The function G{m, x,w^) satisfies the equation 



(mdm - xdx - 2w^d^o'j G{m, x, ~ ~ ^^j G{m, x, = ^ ^ (-^^op - 
Inserting (|B.7| ) in ( |6.6| ) we obtain 

/ j^^dmHQ'^ m)A{,jyp {q; m)r^p(^^)^s (^)^;3 (q) A^s^,^ {q; m) 

= li&^r (^) (^) ^ ^P-/5(np) / (3 + s - 2nz2") z«e-^'"G (m, x, f ) . 

+ 

Considering the leading term [G (m, x, |) = G' (m, x) + O (|)) we have 

= il?- (^)' r (^) r (i±|^) (^) ^ {m, X) + .... 

The Euler's functions T (^ ^+^^~'^ ^ and T ^"^^ °^ order of unity. With the same method 

we prove that (^3) is order O . The previous arguments show a decreasing with the powers of 
the spectrum size of the noise function h{x; m) for the considered loops ordering. This argument 
gives a meaning to the last term of Tab:2 in Section 6. Finally we observe that the previous 
discussion is not a proof, but only a consistency argument with the hypothesis of Section 5. 
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